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1 Introduction 



The origin and evolution of large scale structure is one of the most important problems 
in cosmology today. It is widely accepted that the growth of small fluctuations by 
gravitational instability leads to structure formation. The inflationary paradigm, 
whilst providing possible solutions to a number of other problems associated with the 
hot big bang model, also produces a Gaussian, adiabatic fluctuation spectrum which 
is nearly, though not exactly, scale-invariant (Guth 1981; Albrecht & Steinhardt 1982; 
Linde 1982; Olive 1990; Liddle & Lyth 1993). Based on this prediction, the standard 
Cold Dark Matter (CDM) model of galaxy formation employs the flat, Harrison- 
Zel'dovich spectrum as an input parameter (Efstathiou 1990). The CDM model 
successfully accounts for small (< lOh^^ Mpc) and intermediate {10h~^ Mpc - lOO/i^^ 
Mpc) scale observations, if one introduces a bias in the distribution of luminous to 
dark matter (Davies et al. 1985). Q 

However, standard CDM has come under severe pressure from a number of re- 
cent observations (for a detailed review see Liddle & Lyth 1993). In particular, the 
APM angular galaxy-galaxy correlation function (Maddox et al. 1991) and the IRAS 
QDOT redshift survey (Efstathiou et al. 1991) indicate that there exists more large 
scale structure than that predicted by standard CDM. One possible resolution to 
this problem is to consider tilted CDM models. Here the primordial power spec- 
trum is assumed to be of the form P{k) oc A'g{k)k oc A;", where k is the comoving 
wavenumber of the Fourier expansion of the perturbation, is the amplitude of 
the quantum fluctuation when it crosses the Hubble radius during the matter- or 
radiation-dominated eras and n is the power spectrum. Other possibilities involve 
the addition of a cosmological constant or a hot dark matter component (Liddle & 
Lyth 1993). 

Inflation also produces a spectrum of gravitational wave (tensor) perturbations, 
whose amplitude may or may not be comparable to that of the scalar fluctuations. 
In this paper we shall concentrate on models which lead to tilted power spectra with 
a negligible gravitational wave component. There exists a wide range of observa- 
tional constraints on the tilt arising from large angle {6 > 3°) microwave background 
anisotropics (Smoot et al. 1992), galaxy clustering (Maddox et al. 1990; Efstathiou 
et al. 1990), peculiar velocity flows (Bertschinger & Dekel 1989; Dekel, Bertschinger 
& Faber 1990; Bertschinger et al. 1990), high redshift quasars (Efstathiou & Rees 
1988) and the red shift of structure formation (Adams et al. 1993). When combined 
together these observations strongly limit the allowed value of n. It has been shown 
that tilted CDM can not fit all of the current data simultaneously (Adams et al. 1993; 
Liddle & Lyth 1993). For inflationary models in which gravitational wave production 
is negligible, a lower limit of n > 0.7 is partially consistent with the COBE 2-sigma 
upper limit and the bulk flow data, if the clustering and pairwise velocity data are 
ignored. In models where the gravitational wave contribution to the microwave back- 

^The current value of the expansion rate is Hq = lOOh km s^^ Mpc^^, where 0.4 < h < I. 
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ground anisotropy is important, however, this hmit is strengthened to n > 0.84. This 
is clearly inconsistent with the APM galaxy correlation function, which indicates that 
0.3 < n < 0.7 provides a good fit to the excess clustering data. 

On the other hand, there is growing observational evidence for a departure from 
a pure power law at a scale A ^ 150 ± bOh^^ Mpc (Einasto et al. 1993). The power 
spectrum of clusters of galaxies has spectral index — 2 < n < —1 on intermediate 
scales, whilst it is consistent with n = 1 on large scales (Peacock 1991). 

In short, the current status of the observations is far from conclusive and it is 
therefore important to consider all the theoretical options available. In this paper we 
investigate the general circumstances in which a tilted power law scalar spectrum and 
a negligible gravitational wave amplitude arise in inflationary models. In Section 2, 
we summarize the details of a powerful framework which allows the general form of 
the inflaton potential to be derived in a straightforward manner. The form of such a 
potential is shown to be a hyperbolic secant function in Section 3. In Section 4, it is 
further shown that such a potential arises when a bulk viscous stress is added to the 
energy-momentum tensor of a perfect baryotropic fluid. It is illustrated how a number 
of plausible particle physics models, such as the quantum creation of fundamental 
strings (Turok 1988) and N — 2 supergravity (Salam & Sezgin 1984), lead to a 
potential of this form. 



2 Inflaton dynamics and primordial fluctuations 



Inflation proceeds if the potential energy of a scalar field dominates the energy- 
momentum tensor of the universe, since this leads to a violation of the strong energy 
condition (Hawking & Ellis 1973). One can model self-interacting scalar fields in 
terms of perfect fluids, and vice- versa. Within the context of the isotropic and homo- 
geneous Priedmann-Robcrtson- Walker (FRW) cosmologies, the complete history of 
the universe is determined by three independent equations. The Friedmann equation 

is the first integral of the Raychaudhuri equation (Raychaudhuri 1955) and describes 
the conservation of energy. Here, p(t) represents the total energy density of the 
universe, A; = { — 1, 0,-1-1} for open, flat, and closed spatial sections respectively, a{t) is 
the scale factor, and H = a/ais the expansion rate, where a dot denotes differentiation 
with respect to cosmic time t. Units are chosen such that h = c = SuG = 1. For 
a matter content with pressure p{t), the Bianchi identity (equivalent to the local 
conservation of energy-momentum) is 

p + 3H{p + p)^0, (2) 
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and, in principle, the general solution to these equations follows once the equation of 
state 



P = pip) (3) 

is specified. This equation describes the particle physics sector of the model. 

However, analytical solutions have only been found for a limited number of specific 
examples, such as the special case p/p = constant (Barrow 1990). For an arbitrary 
equation of state which satisfies the dominant energy condition (p — p > 0), it proves 
convenient to redefine the sum and difference of p and p in terms of the new functions 



^'' = p + p ^ (j){t) = J dt^pit') + p{t') (4) 
2V = p-p. (5) 



We may then rewrite equations (|1]) and (|^) as 



which are the Einstein field equations for a minimally coupled scalar field, 0. Specify- 
ing the equation of state now becomes a question of choosing an appropriate functional 
form for the potential V{(j)), and vice- versa. For example, an exponential potential is 
equivalent to p/p = constant when k = 0. 

Recent advances in the treatment of equations (|^) and (|^) have been made by 
viewing the scalar field as the effective dynamical variable of the system (Muslimov 
1990; Salopek & Bond 1990, 1991; Lidsey 1991, 1993). From the definition p = 
(f)'^/2 + V, the time dependence can be eliminated by rewriting the scalar field equation 
(0) as 

p' = -3H^, 0^0, (8) 

where a prime denotes differentiation with respect to 0. This is consistent if does 
not oscillate (i.e. does not pass through zero). It follows that 6H^ = —p'X'/X, 
where X{(f)) = a^(0), and the Friedmann equation becomes 

p'X' + 2pX = Qk. (9) 

The potential can be found immediately from the expression 

once the forms of p(0) and X{(j)) are known. 
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When k = 0, 3H^ = p, and these field equations take the particularly simple form 

2H'a' = -Ha, 2H' = -<j), (11) 

thereby allowing the general solution, a{(f)), to be expressed in terms of quadratures 
with respect to (j) (Salopck & Bond 1991; Lidsey 1993). The expression for the 
potential reduces to an Hamilton- Jacobi differential equation in H{(f)) of the form 

V{(P) = 3//2(0) _ 2{H'y. (12) 

This framework is useful for determining the primordial power spectrum. During 
inflation short-wavelength quantum fluctuations in the inflaton and graviton flelds 
are redshifted beyond the Hubble radius when the comoving wavenumber of the per- 
turbation satisfies k = aH. Once outside, their amplitude remains frozen until they 
re-enter the Hubble radius during the radiation-dominated or matter-dominated eras. 
The density perturbation (scalar) spectrum is given by (Copeland et al. 1993) 

where the quantities on the right are calculated when the fluctuation flrst crosses the 
Hubble radius. In the uniform Hubble constant gauge the constant m = 4 if the 
perturbation re-enters during the radiation-dominated era, whereas m = 2/5 if the 
perturbation re-enters during matter domination. 

The gravitational wave (tensor) spectrum is calculated in the transverse-traceless 
gauge, where -|- and x denote the two independent polarization states of the metric 
perturbation. The classical amplitude of the fluctuation satisfies the massless Klein- 
Gordon equation (Grishchuk 1974, 1977; White 1992), so the graviton can be viewed 
as a massless, minimally coupled scalar field, which has two degrees of freedom 'V'+,x- 
The spectrum of tensor perturbations is then given by (Abbott & Wise 1984) 

One may define two parameters (Copeland et al. 1993) 
and 

.-^^2^. (16) 

H(j) H ^ ^ 

Up to a numerical factor, e measures the relative contribution of the field's kinetic 
energy to its total energy density and may be referred to as the energy parameter. 




4 



The quantity rj measures the ratio of the field's acceleration relative to the friction 
contribution. We refer to it as the friction parameter. 

These parameters may be recast in an alternative form by defining new variables 

x^H^ {H'f. (17) 

It follows from the identity 2H" = d{H'f/dH that 

e = — , V = ^^- (18) 

Infiation proceeds in the region of parameter space for which e < 1 and the coasting 
solution, or Milne universe, corresponds to e = 1 {i.e. 2y = x). It is interesting to 
note that the friction parameter does not directly determine whether infiation occurs. 
The 'slow-roir approximation is valid when {e, \r]\} <^ 1. 

One may also write the (scale-dependent) spectral indices of the scalar and tensor 
fiuctuations in terms of these two quantities. It is easy to show that 

where * indicates that e and rj should be evaluated when a particular scale first crosses 
the Hubble radius. The fiat Harrison-Zel'dovich spectrum is equivalent to n = 1. 
It follows from the definitions of As and Aq that 



^ = ^. (21) 

As m ^ ^ 

It is often stated that infiation leads to a Harrison-Zel'dovich scalar spectrum with a 
negligible gravitational wave contribution. However such a conclusion follows because 
the slow-roll approximation is assumed a priori. Equation ( pi]) implies that the 
gravity wave amplitude can be comparable to As if e is sufficiently large. There 
exists a table of correspondences which summarizes the four possibilities (Barrow & 
Liddle 1993). 

In the following section we shall employ this framework to derive the form of the 
infiaton potential. 



3 The general potential leading to tilted spectra 



Exact results may be derived by solving equation (|TP|). After substitution of equation 
(p!8|) this becomes the first-order differential equation 

dy / n — 5\ y / 1 — n\ . . 



dx \ 2 J X 
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Scalar 
Spectrum 


Gravitational Waves 
Important 


Gravitational Waves 
Negligible 


Small 


e large 


e small 


Tilt 




?7 small 


Significant 


e large 


e small 


Tilt 


?7 large 


\ri\ large 



Table 1 - The table of correspondences illustrating the connection between tilt and the magni- 
tude of the energy and friction parameters. The description 'large' implies significantly larger than 
zero (but still less than unity) and 'small' implies the parameter is very close to zero. 



The solution for constant n is 



1 /I -n 



y = 7; 



2\3-n 



X 



(23) 



where C is an arbitrary integration constant. The case C < is of interest here 
because it leads to the general solution 



ff(0) = A 



sech 



where 



A 



1 n-l 
2|C|n-3 



'(n- l)(n-3) 



l/(3-n) 



2/(3-n) 



n < 1, 



(24) 



(25) 



Without loss of generality the integration constant has been removed by performing 
a linear translation on the value of 0. The constant |C| determines the energy scale 
at which these processes are occurring (Carr & Lidsey 1993; Lidsey & Tavakol 1993). 
From equation (|12D the potential is 



V{(t)) = A^ [sech{u^ 



l4/(3-n) 



n — 1 
n — 3 



tanh (ujc 



where 



(n- l)(n-3) 



The solution to equation 
H{(j)) oc exp 



when C = is the exponential potential 



2\n-3 



(26) 



(27) 



(28) 
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Eqs. (p^ ) and ( pH]) lead to the same scalar perturbation spectrum and a determination 
of As{k) alone will not lift the degeneracy. What distinguishes the two solutions is the 
relative contribution of the gravitational wave spectrum, as determined by equation 
(|1|). We find 



n 



and 



n 



n 



n 



tanh 



'(n- l)(n-3) 



(29) 



(30) 



for the exponential and hyperbolic cases respectively. If cosmological scales crossed 
the Hubble radius when |0| <^ 1, the amplitude of tensor perturbations is exponen- 
tially suppressed in the latter example. 

Since equation (^4]) is an exact solution, it is valid for all values of 0. In particular, 
for sufficiently small (f) the Taylor expansion 



H{<P) = X 



1 



1 



n 



(31) 



will also lead to a constant spectral index. This implies that a power spectrum 
with n = constant < 1 will arise from any function of H{(f)) which is identical to 
equation (|3lD in the small (p approximation. It is well known that deviations from scale 
invariance, without significant gravitational wave production, are possible whenever 
the potential resembles an inverted harmonic oscillator (Steinhardt & Turner 1984), 
but the above calculation provides further insight. Such a result follows because 
the inverted oscillator resembles the hyperbolic secant function to lowest order in 0. 



Equation (^) is very useful because it directly relates the effective imaginary mass 
of the field to the scalar spectral index. 

These results are summarized pictorially in Figs, la and lb, which are representa- 
tions of the class of solutions (P^D in the x — y plane, {x and y are defined in equation 
(p!8|)). In Fig. la the coasting solution y = x/2 is shown as the dashed line and the 
strong energy condition is violated to the right of this line. The x-axis represents the 
de Sitter solution, H = constant, and the origin is Minkowski space, which itself may 
be viewed as de Sitter space with an infinite radius of curvature (Hawking & Ellis 
1973). The solid lines represent solutions of constant n when C = 0. In these models 
X is a measure of the energy density of the universe and decreases as time increases. 
The trajectories of these constant n universes are indicated by the arrows and they 
all asymptotically approach Minkowski space at t +oo.0 

Figure 1 



^In reality the shape of the potential must change at some point {x, y) to allow for an exit from 
inflation. 
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Fig. lb illustrates the trajectories for finite values of C < and n = 0.7. This 
class of universe begins in a de Sitter phase at t = —oo and evolves towards the C = 
asymptote at t = +oo in such a way that the scalar spectral index remains constant 
at all times. The magnitude of C determines the amplitude of the scalar quantum 
fluctuations but not the scale dependence. 

Having found the form of the potential required, it is now necessary to consider 
the physics which may lead to such a model. 



4 Modelling the potential as a bulk viscous stress 

In this section we employ the techniques discussed in Section 2 for an equation of 
state which has a number of physical applications in the early universe. Equation (^) 
may be rewritten in the form p = T{p)p for some arbitrary function F. The simplest 
form for F(p) is the baryotropic equation of state 

r = 7 - 1 = constant, < 7 < 2, (32) 

which is equivalent at the classical level to the subset of solutions 

0' = ^^V, /X ^ (33) 

2-7 

Hence < yU < 00, where = (7 = 0) corresponds to vacuum energy and /i = 
cxD (7 = 2) represents a stiff fluid. In particular, radiation and matter dominated 
universes are characterized by /i = 4 (7 = 4/3) and /i = 2 (7 = 1) respectively. 
Furthermore, the equation of state for a universe dominated by topological defects 
such as domain walls or a gas of cosmic strings is also given by equation (|32|) . It 
can be shown that 714/ = v'y/ + 1/3 for a gas of domain walls moving with constant 
velocity vw^ whereas 75 = 2(1 + v'g)/?) for a string gas (Kolb & Turner 1990). 

Mathematically, a natural and plausible extension towards a more realistic effec- 
tive equation of state is to view the function F as a truncated Taylor series 

F(p) = (7-l)-/5p°, (34) 

for positive definite constants {a,j3}. The second term in this expression can be 
treated as a first-order perturbative correction to F = constant. Physically, when 
= 0, equation (|3^) is equivalent to a fluid with bulk viscosity ?7(p) = Hp*" and a 
perfect fluid p = (7 — l)p, where (3 = a/3H and m = a + \ (Barrow 1988). It is 
known that the effects of particle production in the early universe can be modelled in 
terms of a classical bulk viscosity of this form. The polarization (trace anomaly) of 
quantized fields in curved spacetimes also leads to vacuum viscosity effects (Zel'dovich 
1980; Hu 1982; Waga, Falcao & Chanda 1986). Exact cosmological solutions based 
on equation (R^ have been found by Murphy (1973) and Barrow (1988), but in this 
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section we derive them in a much simpler parametric form. This allows the primordial 
fluctuation spectra to be calculated. 

With this choice of r(p) the Bianchi identity integrates exactly for all k to 
yield 



p(a) 



7 V«o/ 



-l/a 



(35) 



where the constant of integration is expressed in terms of ao and {0,7} 7^ 0. When 
/3 = 0, p oc a'^"^ as expected. Unfortunately analytical solutions to the Friedmann 
equation (|^) have not been found when k = ±1, but this equation simplifies to 

^ = 3 (7 - /3p") (36) 
P 

when k = and this has the exact solution 

1 - -p°(0) = tanh2(cj0) (37) 
7 

. . ,38) 

The dominant energy condition is violated when /5 < 0, so this case is not considered 
further. The expressions for if(0), a(0) and V{(j)) follow from equations (0) and (|12[) 

as 



W) = ^(^^l [sech(.;0)]V" (39) 



/ n\ 1/307 

a(0) = ao i^j |sinh(cu0)|2/'"^ (40) 



7 ^ 



^(</>) = (^^j [sech{u<j,)f^ (1 - (7/2)tanh2(^</))) . (41) 

It should be emphasized that these solutions are exact and no 'slow-roll' approxima- 
tions, such as 101 <C H\(j)\ and 0^ -C V, have been made. These parametric solutions 
are plotted schematically in Fig. 2. 

Figure 2 

For completeness, we include the well known solutions for /? = 0, which lead to 
the exponential potential 

H{^) (X exp h/ ^0 I , V(0) cx exp ( J370 ) (42) 
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and power law expansion a(t) oc t'^^^"'. 

Hence, when = 0, a perfect baryotropic fluid with bulk viscosity can be modelled 



as a self-interacting scalar field with potential (|T). The inclusion of bulk viscosity 
alters the structure of the potential away from an exponential form in the neighbour- 
hood of |0| ~ 0. This is illustrated in Fig. 2c. Near the origin, we find that the 
equation of state (0) can be adequately described as an inverted harmonic oscillator. 
The exponential potentials are recovered in the asymptotic limit as |0| oo, because 
the viscous effects decay faster than the perfect fluid contribution at large |0|. The 
general advantage of rewriting equation in terms of a scalar field is that the qual- 
itative history of the universe is easily determined by considering the evolution of the 
scalar field along its interaction potential. Here the field is initially placed at = 0, 
which corresponds to a de Sitter expansion with H = (7//3)^/2a/y3_ the field rolls 
either to the left or right of the potential, the universe expands and approaches the 
power law attractor solution a oc t'^^^"' with H monotonically. This behaviour 
was first noted by Barrow (1988). 

The region of parameter space in which the strong energy condition (p + 3p > 0) 
is violated can be determined from equations ( ^41) and (^Tf). Defining the quantity 
A = p + 3p, we find 

A = 37P |^tanh^(tu0) - . (43) 

It follows that inflation occurs for all (p when 7 < 2/3 and there exists a graceful exit 
problem. An identical problem is encountered when (3 = 0. The potential must be 
modifled in such a way that allows the expansion to become sub luminal. For 7 > 2/3, 
on the other hand, the universe inflates initially, but deflates once 

2 

tanh^cj^ > — (44) 
37 

is satisfled. The end of inflation can be deflned as the point where this becomes an 
equality. 

The epoch during inflation which is relevant for large-scale structure observations 
can then be determined and we now proceed to discuss models which lead to a po- 
tential of the form given by equation 



5 Particle physics models 

The purpose of this section is to discuss a number of particle physics models in which 
the features described above may arise. 

A. Bulk Viscosity 
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A direct comparison of Eqs. (p4D and (|39D indicates that tilted spectra with 



constant spectral index and negligible gravitational waves arise in a class of bulk 
viscosity models if we identify 

" = ^ (45) 

' 3\n-3j 37-2' ^ ' 

or equivalently 

The constraint n > 0.7 now becomes the upper limit 7 < 0.09. (We require 7 < 2/3 
for consistency in this class, since equation (|39D is only valid for positive a). The 
parameter f3 of equation (|3^) plays the role of \C\ in determining the amplitude of 
the fluctuations, whereas a determines the tilt of the spectrum. 

When {a, 7} are not related by equation ( ^Tj) the tilt is not exactly scale invariant, 
but n ~ constant is an excellent approximation if Lj\(j)\ -C 1. For standard reheating, 
scales of astrophysical interest first crossed the Hubble radius approximately 50 e- 
foldings before the end of inflation. Defining the value of the field at this point as 
050, we find from equation (|^) that 

smh\uM = (^:^) e"'"^"^" (48) 

for 7 > 2/3, where N^o ~ 50. Hence, if 7 is not too close to 7 = 2/3, ujcp^o <^ 1 
is valid and it is consistent to expand equation (RDI) to lowest order. It follows from 
equation (|3ll) that 

n = l-3aj (49) 
and n > 0.7 leads to the constraint 

07 < 0.1. (50) 



B. Quantum creation of fundamental strings 

Turok (1988) has considered the quantum production of infinitely thin Witten 
strings on super- horizon size scales (Green, Schwarz & Witten 1988). He suggested 
that a deflationary expansion follows naturally from a quasi de Sitter phase in the 
early universe. When a = 1 and 2/3<7<l,it was shown that the equation of state 
(0) provides a good phenomenological description of the quantum creation of these 
strings after compactification to four dimensions (Barrow 1988; Turok 1988). The 
parameter (3 depends on the fundamental string tension and the fractal dimension of 
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the string. The lower and upper hmits on 7 correspond to long strings with negligible 
velocity and a highly convoluted, relativistic string distribution respectively. 

Therefore, the parametric solutions (|39|)-PTD describe the evolution of the flat 
Friedmann universe when dominated by fundamental strings created on super-horizon 
scales. It follows from equation ( P9| ) that —2 < n < — 1 and the model in its present 
form gives far too much power on large scales. However, it is interesting that this 
predicted range for the spectral index corresponds precisely to that observed on inter- 
mediate scales 1 < A < 200/i^^ Mpc in the distribution of clusters of galaxies (Einasto 
et al. 1993). 

C. = 2 supergravity in six dimensions 

The Salam-Sezgin model is = 2 supergravity in six dimensions compactified 
onto a two-sphere (Salam & Sezgin 1984). Although such a model does not reproduce 
a correct particle spectrum in four dimensions, it is thought to contain features generic 
to more complete theories and its cosmological implications have been studied in some 
detail (Liddle 1989). It has been shown that the theory is equivalent to two interacting 
scalar fields, ^ and a, whose full potential exhibits a global minimum in both ^ and a 
directions if the scale-invariance of the theory is broken by quantum effects (Gibbons 
& Townsend 1987). If the ^-field comes to rest, the potential for a becomes 

V{a) = Vo exp(— v^cr) — 2exp i-V2a) + 1 . (51) 

The potential has a global minimum located at a = and an asymptotic form V cx 
exp(— a/Scj) for a <^ 0. In this regime the system behaves as a baryotropic perfect 
fluid with 7 = 8/3. In this unmodified form the model does not lead to inflation, but 
it is clear from Sections 3 and 4 that introducing a bulk viscosity into the system will 
violate the strong energy condition and lead to a tilted spectrum. By condition ( ^OD 
observation requires a < 0.038. 

Moreover, the viscosity effects redshift faster than the perfect fluid contribution, 
implying that they will become negligible by the time the field has rolled towards 
0" = 0. Hence, the shape of the potential is not significantly altered within the vicinity 
of the global minimum and standard reheating can occur through rapid oscillations 
of the field about this minimum (Kolb & Turner 1990). This solves the graceful exit 
problem discussed in Section 2. 

D. Natural Inflation 

For completeness we remark that natural inflation, driven by a pseudo-Nambu- 
Goldstone boson with potential V{(f)) oc H-cos(0//), also produces a tilted spectrum. 
This model has been discussed in detail by Adams et al. (1993). In the small angle 
approximation the spectral index is given by n = 1 — via equation (pl|). 
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6 Conclusions 



The general inflationary potential which produces a measureable tilt in the primor- 
dial spectrum of scalar fluctuations, without producing a significant gravitational 
wave spectrum, is of the form (|26|). In this paper it has been shown that such a po- 
tential is physically equivalent to a bulk viscous stress modification to a baryotropic 
perfect fluid. Mathematically it is equivalent to expanding the equation of state as a 
truncated Taylor series. It was further shown that this potential arises in a number of 
particle physics theories. If future observations indicate that a) gravitational waves 
are not contributing to the anisotropy in the cosmic microwave background and b) a 
dark matter model based on a spectrum with significant tilt does provide a good fit 
to the data, then this would justify a more detailed study of the models discussed in 
Section 5. 

A number of simplifying assumptions were made. In particular, the expressions 
(p!3|) and (|1^ for the scalar and tensor amplitudes are strictly only valid in the slow 
roll regime, {e, \ri\} <^ 1, whereas equation (|T^) implies that a tilt away from the 
Harrison-Zel'dovich spectrum requires <^ Ir^l < 1, approximately. However, it has 
been shown that the corrections away from slow-roll are not important near a local 
maximum, and although they slightly alter the amplitude of the fluctuations in the 
exponential regime, they do not change the spectral index (Stewart & Lyth 1993). 

The physical interpretation of the potential in terms of a bulk viscosity is only 
valid in the spatially flat FRW cosmology. However, this paper has investigated the 
power spectra of such models and it is the last 60 e-foldings of inflationary expansion 
which are important for large-scale structure (Kolb & Turner 1990). In most chaotic 
scenarios the density parameter is very close to unity by this stage. 

Furthermore, since all scales probed by large-scale structure correspond to a small 
(~ 9) number of e-foldings, it is reasonable to assume that the parameter 7 is constant 
during this interval. These results may therefore have wider applications in models 
where the polytropic index is a function of cosmic time. 

Finally we note that equation with 7 = 1 and /5 < is the equation of state for 
a polytropic star, special cases of which include white dwarfs (a = 5/3) and neutron 
stars (a = 4/3) (Weinberg 1973). If the solution to equation (|^) could be found for 
k = +1, the techniques described here could be relevant for stellar structure. 
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Figure Captions 



Figure la: The trajectories of the constant n universes in the x — y place when 
C = 0. The dashed hne represents the Milne universe and the strong energy condition 
is satisfied to the left of this line. All trajectories approach the origin at t = +00, 
which is de Sitter space with infinite radius of curvature. 

Figure lb: The trajectories of n = 0.7 universes when C = {0, —0.5, —1, —2}. The 
dashed line corresponds to C = 0. All solutions begin in a de Sitter phase at t = —00 
and approach the C — (power law) asymptote at t — +00. 

Figure 2: Schematic plots of the parametric solution (38)- (40) as a function of 
the scalar field. In Fig. 2c the dashed lines represent the exponential potential corre- 
sponding to vanishing bulk viscosity. The inclusion of a bulk viscous stress results in 
a local maximum in the potential and a quasi de Sitter expansion. It is this feature 
which results in negligible gravitational wave production, but the curvature of the 
potential at the origin tilts the scalar spectrum. 
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